DeepSub: Utilizing Deep Learning for Predicting the Number of Subunits in Homo-Oligomeric Protein Complexes.

The molecular weight (MW) of an enzyme is a critical parameter in enzyme-constrained models (ecModels). It is determined by two factors: the presence of subunits and the abundance of each subunit. Although the number of subunits (NS) can potentially be obtained from UniProt, this information is not readily available for most proteins. In this study, we addressed this gap by extracting and curating subunit information from the UniProt database to establish a robust benchmark dataset. Subsequently, we propose a novel model named DeepSub, which leverages the protein language model and Bi-directional Gated Recurrent Unit (GRU), to predict NS in homo-oligomers solely based on protein sequences. DeepSub demonstrates remarkable accuracy, achieving an accuracy rate as high as 0.967, surpassing the performance of QUEEN. To validate the effectiveness of DeepSub, we performed predictions for protein homo-oligomers that have been reported in the literature but are not documented in the UniProt database. Examples include homoserine dehydrogenase from Corynebacterium glutamicum, Matrilin-4 from Mus musculus and Homo sapiens, and the Multimerins protein family from M. musculus and H. sapiens. The predicted results align closely with the reported findings in the literature, underscoring the reliability and utility of DeepSub.

Analysis of Surface Expression of NMDAR Subunits in Primary Hippocampal Neurons.

The expression and activity of ionotropic glutamate receptors control signal transduction at the excitatory synapses in the CNS. The NMDAR comprises two obligatory GluN1 subunits and two GluN2 or GluN3 subunits in different combinations. Each GluN subunit consists of four domains: the extracellular amino-terminal and agonist-binding domains, the transmembrane domain, and the intracellular C-terminal domain (CTD). The CTD interaction with various classes of intracellular proteins is critical for trafficking and synaptic localization of NMDARs. Amino acid mutations or the inclusion of premature stop codons in the CTD could contribute to the emergence of neurodevelopmental and neuropsychiatric disorders. Here, we describe the method of preparing primary hippocampal neurons and lentiviral particles expressing GluN subunits that can be used as a model to study cell surface expression and synaptic localization of NMDARs. We also show a simple method of fluorescence immunostaining of eGFP-tagged GluN2 subunits and subsequent microscopy technique and image analysis to study the effects of disease-associated mutations in the CTDs of GluN2A and GluN2B subunits.

Functional Analysis of NMDAR Subunit Components in Postsynaptic Currents of Identified Cells and Synapses in Brain Slices.

Epilepsy is one of the most represented neurological diseases worldwide. However, in many cases, the precise molecular mechanisms of epileptogenesis and ictiogenesis are unknown. Because of their important role in synaptic function and neuronal excitability, NMDA receptors are implicated in various epileptogenic mechanisms. Most of these are subunit specific and require a precise analysis of the subunit composition of the NMDARs implicated. Here, we describe an express electrophysiological method to analyze the contribution of NMDAR subunits to spontaneous postsynaptic activity in identified cells in brain slices using patch clamp whole cell recordings.

Selective Cell-Surface Expression of Triheteromeric NMDA Receptors.

NMDA-type ionotropic glutamate receptors are critically involved in many brain functions and are implicated in a variety of brain disorders. Seven NMDA receptor subunits exist (GluN1, GluN2A-D, and GluN3A-B) that assemble into tetrameric receptor subtypes with distinct functional properties and physiological roles. The majority NMDA receptors are composed of two GluN1 and two GluN2 subunits, which can assemble into four diheteromeric receptors subtypes composed of GluN1 and one type of GluN2 subunit (e.g., GluN1/2A), and presumably also six triheteromeric receptor subtypes composed of GluN1 and two different GluN2 subunits (e.g., GluN1/2A/2B). Furthermore, the GluN1 subunit exists as eight splice variants (e.g., GluN1-1a and GluN1-1b isoforms), and two different GluN1 isoforms can co-assemble to also form triheteromeric NMDA receptors (e.g., GluN1-1a/1b/2A). Here, we describe a method to faithfully express triheteromeric NMDA receptors in heterologous expression systems by controlling the identity of two of the four subunits. This method overcomes the problem that co-expression of three different NMDA receptor subunits generates two distinct diheteromeric receptor subtypes as well as one triheteromeric receptor subtype, thereby confounding studies that require a homogenous population of triheteromeric NMDA receptors. The method has been applied to selectively express recombinant triheteromeric GluN1/2A/2B, GluN1/2A/2C, GluN1/2B/2D, GluN1-1a/GluN1-1b/2A, GluN1-1a/GluN1-1b/2B receptors with negligible co-expression of the respective diheteromeric receptor subtypes. This method therefore enables quantitative evaluation of functional and pharmacological properties of triheteromeric NMDA receptors, some of which are abundant NMDA receptor subtypes in the adult brain.

Characterizing Functional Contributions of Specific GluN2 Subunits to Individual Postsynaptic NMDAR Responses Using Biophysical Parameters.

The NMDAR is a heterotetramer composed of two GluN1 subunits and two GluN2 and/or GluN3 subunits, with the GluN2 subunits exhibiting significant diversity in their structure and function. Recent studies have highlighted the importance of characterizing the specific roles of each GluN2 subunit across central nervous system regions and developmental stages, as well as their unique contributions to NMDAR-mediated signaling and plasticity. Understanding the distinct functions of GluN2 subunits is critical for the development of targeted therapeutic strategies for NMDAR-related disorders. However, measuring the functional contribution of individual GluN2 subtypes in ex vivo slices is challenging. Conventionally, pharmacological or genetic approaches are used, but, in many cases, this is not possible or is restricted to population-level NMDAR responses. Here, we describe a technique for using biophysical properties of miniature synaptic NMDAR responses as a proxy to measure the functional contribution of specific GluN2-NMDAR subunits to individual synapses within a neuron.

CORVET-specific subunit levels determine the balance between HOPS/CORVET endosomal tethering complexes.

The closely related endolysosomal tethering complexes HOPS and CORVET play pivotal roles in the homo- and heterotypic fusion of early and late endosomes, respectively, and HOPS also mediates the fusion of lysosomes with incoming vesicles including late endosomes and autophagosomes. These heterohexameric complexes share their four core subunits that assemble with additional two, complex-specific subunits. These features and the similar structure of the complexes could allow the formation of hybrid complexes, and the complex specific subunits may compete for binding to the core. Indeed, our biochemical analyses revealed the overlap of binding sites for HOPS-specific VPS41 and CORVET-specific VPS8 on the shared core subunit VPS18. We found that the overexpression of CORVET-specific VPS8 or Tgfbrap1 decreased the amount of core proteins VPS11 and VPS18 that are assembled with HOPS-specific subunits VPS41 or VPS39, indicating reduced amount of assembled HOPS. In line with this, we observed the elevation of both lipidated, autophagosome-associated LC3 protein and the autophagic cargo p62 in these cells, suggesting impaired autophagosome-lysosome fusion. In contrast, overexpression of HOPS-specific VPS39 or VPS41 did not affect the level of assembled CORVET or autophagy. VPS8 or Tgfbrap1 overexpression also induced Cathepsin D accumulation, suggesting that HOPS-dependent biosynthetic delivery of lysosomal hydrolases is perturbed, too. These indicate that CORVET-specific subunit levels fine-tune HOPS assembly and activity in vivo.

Cryo-EM Structure of Porphyromonas gingivalis RNA Polymerase.

Porphyromonas gingivalis, an anaerobic CFB (Cytophaga, Fusobacterium, and Bacteroides) group bacterium, is the keystone pathogen of periodontitis and has been implicated in various systemic diseases. Increased antibiotic resistance and lack of effective antibiotics necessitate a search for new intervention strategies. Here we report a 3.5 Å resolution cryo-EM structure of P. gingivalis RNA polymerase (RNAP). The structure displays new structural features in its ω subunit and multiple domains in ß and ß' subunits, which differ from their counterparts in other bacterial RNAPs. Superimpositions with E. coli RNAP holoenzyme and initiation complex further suggest that its ω subunit may contact the σ4 domain, thereby possibly contributing to the assembly and stabilization of initiation complexes. In addition to revealing the unique features of P. gingivalis RNAP, our work offers a framework for future studies of transcription regulation in this important pathogen, as well as for structure-based drug development.

Alternative splicing for tuneable expression of protein subunits at desired ratios.

The controlled expression of two or more proteins at a defined and stable ratio remains a substantial challenge, particularly in the bi- and multispecific antibody field. Achieving an optimal ratio of protein subunits can facilitate the assembly of multimeric proteins with high efficiency and minimize the production of by-products. In this study, we propose a solution based on alternative splicing, enabling the expression of a tunable and predefined ratio of two distinct polypeptide chains from the same pre-mRNA under the control of a single promoter. The pre-mRNA used in this study contains two open reading frames situated on separate exons. The first exon is flanked by two copies of the chicken troponin intron 4 (cTNT-I4) and is susceptible to excision from the pre-mRNA by means of alternative splicing. This specific design enables the modulation of the splice ratio by adjusting the strength of the splice acceptor. To illustrate this approach, we developed constructs expressing varying ratios of GFP and dsRED and extended their application to multimeric proteins such as monoclonal antibodies, achieving industrially relevant expression levels (>1 g/L) in a 14-day fed-batch process. The stability of the splice ratio was confirmed by droplet digital PCR in a stable pool cultivated over a 28-day period, while product quality was assessed via intact mass analysis, demonstrating absence of product-related impurities resulting from undesired splice events. Furthermore, we showcased the versatility of the construct by expressing two subunits of a bispecific antibody of the BEAT® type, which contains three distinct subunits in total.

A Method to Conditionally Measure Target Engagement at Intracellular RAS and RAF Complexes.

Dysfunction of the RAS/mitogen-activated protein kinase (MAPK) pathway is a common driver of human cancers. As such, both the master regulator of the pathway, RAS, and its proximal kinase effectors, RAFs, have been of interest as drug targets for decades. Importantly, signaling within the RAS/MAPK pathway is highly coordinated due to the formation of a higher-order complex called the RAS/RAF signalosome, which may minimally contain dimers of both RAS and RAF protomers. In the disease state, RAS and RAF assemble in homo- and/or heterodimeric forms. Traditionally, drug development campaigns for both RAS and RAF have utilized biochemical assays of purified recombinant protein. As these assays do not query the RAS or RAF proteins in their full-length and complexed forms in cells, potency results collected using these assays have often failed to correlate with inhibition of the MAPK pathway. To more accurately quantify engagement at this signaling components, we present a bioluminescence resonance energy transfer (BRET)-based method to conditionally measure target engagement at individual protomers within the RAS/RAF signalosome in live cells.

How an intrinsically disordered regulatory subunit assembles a PP1:eIF2 complex.

Fatalska et al.1 use an interdisciplinary strategy to elucidate how an intrinsically disordered regulatory subunit of protein phosphatase 1 binds trimeric eIF2 and positions the phosphatase-substrate complex for dephosphorylation. As validation, they show that a disease mutation abolishes the interaction.

Structural basis of Integrator-dependent RNA polymerase II termination.

The Integrator complex can terminate RNA polymerase II (Pol II) in the promoter-proximal region of genes. Previous work has shed light on how Integrator binds to the paused elongation complex consisting of Pol II, the DRB sensitivity-inducing factor (DSIF) and the negative elongation factor (NELF) and how it cleaves the nascent RNA transcript1, but has not explained how Integrator removes Pol II from the DNA template. Here we present three cryo-electron microscopy structures of the complete Integrator-PP2A complex in different functional states. The structure of the pre-termination complex reveals a previously unresolved, scorpion-tail-shaped INTS10-INTS13-INTS14-INTS15 module that may use its 'sting' to open the DSIF DNA clamp and facilitate termination. The structure of the post-termination complex shows that the previously unresolved subunit INTS3 and associated sensor of single-stranded DNA complex (SOSS) factors prevent Pol II rebinding to Integrator after termination. The structure of the free Integrator-PP2A complex in an inactive closed conformation2 reveals that INTS6 blocks the PP2A phosphatase active site. These results lead to a model for how Integrator terminates Pol II transcription in three steps that involve major rearrangements.

1H, 13C, and 15N backbone and methyl group resonance assignments of ricin toxin A subunit.

Ricin is a potent plant toxin that targets the eukaryotic ribosome by depurinating an adenine from the sarcin-ricin loop (SRL), a highly conserved stem-loop of the rRNA. As a category-B agent for bioterrorism it is a prime target for therapeutic intervention with antibodies and enzyme blocking inhibitors since no effective therapy exists for ricin. Ricin toxin A subunit (RTA) depurinates the SRL by binding to the P-stalk proteins at a remote site. Stimulation of the N-glycosidase activity of RTA by the P-stalk proteins has been studied extensively by biochemical methods and by X-ray crystallography. The current understanding of RTA's depurination mechanism relies exclusively on X-ray structures of the enzyme in the free state and complexed with transition state analogues. To date we have sparse evidence of conformational dynamics and allosteric regulation of RTA activity that can be exploited in the rational design of inhibitors. Thus, our primary goal here is to apply solution NMR techniques to probe the residue specific structural and dynamic coupling active in RTA as a prerequisite to understand the functional implications of an allosteric network. In this report we present de novo sequence specific amide and sidechain methyl chemical shift assignments of the 267 residue RTA in the free state and in complex with an 11-residue peptide (P11) representing the identical C-terminal sequence of the ribosomal P-stalk proteins. These assignments will facilitate future studies detailing the propagation of binding induced conformational changes in RTA complexed with inhibitors, antibodies, and biologically relevant targets.

Architecture and activation of human muscle phosphorylase kinase.

The study of phosphorylase kinase (PhK)-regulated glycogen metabolism has contributed to the fundamental understanding of protein phosphorylation; however, the molecular mechanism of PhK remains poorly understood. Here we present the high-resolution cryo-electron microscopy structures of human muscle PhK. The 1.3-megadalton PhK α4ß4γ4δ4 hexadecamer consists of a tetramer of tetramer, wherein four αßγδ modules are connected by the central ß4 scaffold. The α- and ß-subunits possess glucoamylase-like domains, but exhibit no detectable enzyme activities. The α-subunit serves as a bridge between the ß-subunit and the γδ subcomplex, and facilitates the γ-subunit to adopt an autoinhibited state. Ca2+-free calmodulin (δ-subunit) binds to the γ-subunit in a compact conformation. Upon binding of Ca2+, a conformational change occurs, allowing for the de-inhibition of the γ-subunit through a spring-loaded mechanism. We also reveal an ADP-binding pocket in the ß-subunit, which plays a role in allosterically enhancing PhK activity. These results provide molecular insights of this important kinase complex.

Structure and assembly of a bacterial gasdermin pore.

In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1-3. Studies of human and mouse GSDM pores have revealed the functions and architectures of assemblies comprising 24 to 33 protomers4-9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing more than 50 protomers. We determine a cryo-electron microscopy structure of a Vitiosangium bGSDM in an active 'slinky'-like oligomeric conformation and analyse bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, our results support a stepwise model of GSDM pore assembly and suggest that a covalently bound palmitoyl can leave a hydrophobic sheath and insert into the membrane before formation of the membrane-spanning ß-strand regions. These results reveal the diversity of GSDM pores found in nature and explain the function of an ancient post-translational modification in enabling programmed host cell death.

RrA, an enzyme from Rhodospirillum rubrum, is a prototype of a new family of short-chain L-asparaginases.

L-Asparaginases (ASNases) catalyze the hydrolysis of L-Asn to L-Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300-400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum, reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N-terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L-Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short-chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild-type RrA, both with and without bound L-Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand-alone subfamily of ASNases that can catalyze the hydrolysis of L-Asn to L-Asp despite the lack of the C-terminal domain that is present in all ASNases described structurally to date.

Exploring the structural assembly of rice ADP-glucose pyrophosphorylase subunits using MD simulation.

ADP-glucose pyrophosphorylase plays a pivotal role as an allosteric enzyme, essential for starch biosynthesis in plants. The higher plant AGPase comparises of a pair of large and a pair of small subunits to form a heterotetrameric complex. Growing evidence indicates that each subunit plays a distinct role in regulating the underlying mechanism of starch biosynthesis. In the rice genome, there are four large subunit genes (OsL1-L4) and three small subunit genes (OsS1, OsS2a, and OsS2b). While the structural assembly of cytosolic rice AGPase subunits (OsL2:OsS2b) has been elucidated, there is currently no such documented research available for plastidial rice AGPases (OsL1:OsS1). In this study, we employed protein modeling and MD simulation approaches to gain insights into the structural association of plastidial rice AGPase subunits. Our results demonstrate that the heterotetrameric association of OsL1:OsS1 is very similar to that of cytosolic OsL2:OsS2b and potato AGPase heterotetramer (StLS:StSS). Moreover, the yeast-two-hybrid results on OsL1:OsS1, which resemble StLS:StSS, suggest a differential protein assembly for OsL2:OsS2b. Thus, the regulatory and catalytic mechanisms for plastidial AGPases (OsL1:OsS1) could be different in rice culm and developing endosperm compared to those of OsL2:OsS2b, which are predominantly found in rice endosperm.

Cryo-EM structures of RAD51 assembled on nucleosomes containing a DSB site.

RAD51 is the central eukaryotic recombinase required for meiotic recombination and mitotic repair of double-strand DNA breaks (DSBs)1,2. However, the mechanism by which RAD51 functions at DSB sites in chromatin has remained elusive. Here we report the cryo-electron microscopy structures of human RAD51-nucleosome complexes, in which RAD51 forms ring and filament conformations. In the ring forms, the N-terminal lobe domains (NLDs) of RAD51 protomers are aligned on the outside of the RAD51 ring, and directly bind to the nucleosomal DNA. The nucleosomal linker DNA that contains the DSB site is recognized by the L1 and L2 loops-active centres that face the central hole of the RAD51 ring. In the filament form, the nucleosomal DNA is peeled by the RAD51 filament extension, and the NLDs of RAD51 protomers proximal to the nucleosome bind to the remaining nucleosomal DNA and histones. Mutations that affect nucleosome-binding residues of the RAD51 NLD decrease nucleosome binding, but barely affect DNA binding in vitro. Consistently, yeast Rad51 mutants with the corresponding mutations are substantially defective in DNA repair in vivo. These results reveal an unexpected function of the RAD51 NLD, and explain the mechanism by which RAD51 associates with nucleosomes, recognizes DSBs and forms the active filament in chromatin.

Assembly mechanisms of the neuronal gap junction channel connexin 36 elucidated by Cryo-EM.

Electrical synapses are essential components of neural circuits. Neuronal signal transduction across electrical synapses is primarily mediated by gap junction channels composed of Connexin36 (Cx36), the lack of which causes impaired electrical coupling between certain neurons including cortical interneurons and thalamic reticular nucleus (TRN) neurons. However, the structural basis underlying Cx36 function and assembly remains elusive. Recently, Lee et al. reported cryo-EM structures of Cx36, thus provided first insights of its gating mechanism. Here, we report a consistent cryo-EM structure of Cx36 determined in parallel, and describe unique interactions underpinning its assembly mechanism in complementary to the competing work. In particular, we found non-canonical electrostatic interactions between protomers from opposing hemichannels and a steric complementary site between adjacent protomers within a hemichannel, which together provide a structural explanation for the assembly specificity in homomeric and heteromeric gap junction channels.

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations.

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Identification of Unique Fragmentation Patterns of Fentanyl Analog Protomers Using Structures for Lossless Ion Manipulations Ion Mobility-Orbitrap Mass Spectrometry.

The opioid crisis in the United States is being fueled by the rapid emergence of new fentanyl analogs and precursors that can elude traditional library-based screening methods, which require data from known reference compounds. Since reference compounds are unavailable for new fentanyl analogs, we examined if fentanyls (fentanyl + fentanyl analogs) could be identified in a reference-free manner using a combination of electrospray ionization (ESI), high-resolution ion mobility (IM) spectrometry, high-resolution mass spectrometry (MS), and higher-energy collision-induced dissociation (MS/MS). We analyzed a mixture containing nine fentanyls and W-15 (a structurally similar molecule) and found that the protonated forms of all fentanyls exhibited two baseline-separated IM distributions that produced different MS/MS patterns. Upon fragmentation, both IM distributions of all fentanyls produced two high intensity fragments, resulting from amine site cleavages. The higher mobility distributions of all fentanyls also produced several low intensity fragments, but surprisingly, these same fragments exhibited much greater intensities in the lower mobility distributions. This observation demonstrates that many fragments of fentanyls predominantly originate from one of two different gas-phase structures (suggestive of protomers). Furthermore, increasing the water concentration in the ESI solution increased the intensity of the lower mobility distribution relative to the higher mobility distribution, which further supports that fentanyls exist as two gas-phase protomers. Our observations on the IM and MS/MS properties of fentanyls can be exploited to positively differentiate fentanyls from other compounds without requiring reference libraries and will hopefully assist first responders and law enforcement in combating new and emerging fentanyls.